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Glueballs are particles whose valence degrees of freedom are gluons and therefore in their descrip- 
tion the gauge field plays a dominant role. We review recent results in the physics of glueballs with 
the aim set on phenomenology and discuss the possibility of finding them in conventional hadronic 
experiments and in the Quark Gluon Plasma. In order to describe their properties we resort to a va- 
riety of theoretical treatments which include, lattice QCD, constituent models, AdS/QCD methods, 
and QCD sum rules. The review is supposed to be an informed guide to the literature. Therefore, 
we do not discuss in detail technical developments but refer the reader to the appropriate references. 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) is the theory of the hadronic interactions. It is an elegant theory whose full 
non perturbative solution has escaped our knowledge since its formulation more than 30 years ago.[l| The theory is 
asymptotically freed HI and confining. [1] A particularly good test of our understanding of the nonperturbative aspects 
of QCD is to study particles where the gauge field plays a more important dynamical role than in the standard hadrons. 
In particular glueballs, bound states of gluons, represent such a scenario. 

The glueball spectrum has attracted much attention since the early days of QCD.p] The interest in this subject 
is related to the significant progress in the understanding of the properties of such states within QCD, as well as, 
in the new possibilities for their identification in modern experiments. Historically the investigation of the glueball 
properties started in the bag model by Jaffe and Johnson. [6] They found many glueball states with different quantum 
numbers lying in the mass interval 1000-2000 MeV. They emphasized that one should expect rather small widths for 
such states because their decays in conventional hadrons violate the Okubo-Zweig-Iizuka (OZI) rule.Q 

After this pioneering work the study of glueballs was carried out by using various versions of constituent models, by 
exploiting the QCD sum rule approach and by performing lattice QCD calculations. Glueballs have not been an easy 
subject to study due to the lack of phcnomcnological support and therefore much debate has been associated with 
their properties. The main achievement of these approaches is the understanding of the deep relation between the 
properties of the glueball states and the structure of the QCD vacuum. Besides, they provide a determination of the 
spectrum both in gluodynamics, the theory with just gluons and no quarks, and in QCD. However, in (unquenched) 
QCD, the results of several calculations for the spectrum are still not universally accepted, in particular, for the lowest 
lying glueballs. US] 

From the phenomenological point of view it has become clear by now that it is difficult to single out which states 
of the hadronic spectrum are glueballs because we lack the necessary knowledge to determine their decay properties. 
Moreover the strong expected mixing between glueballs and quark states leads to a broadening of the possible glueball 
states which does not simplify their isolation. The wishful sharp resonances which would confer the glueball spectra 
the beauty and richness of the baryonic and mesonic spectra are lacking. This confusing picture has led to a loss of 
theoretical and experimental interest in these hadronic states. However, it is important to stress, that if they were to 
exist they would be a beautiful and unique consequence of QCD. At the present, several candidates for the low mass 
glueballs with quantum numbers ++ , 2 ++ , ^ and 1" are under discussion. [Tol. ITTI. Il2l. [T3L [Til HH 
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In this review we will discuss the modern development in glueball spectroscopy from various perspectives. In 
section |n] we will summarize the results that lattice techniques have obtained for the spectrum, both in the pure gauge 
theory and in the unquenched calculations. In section ITlTl we present a review on constituent models. Section Hvl is 
dedicated to discuss QCD sum rules. In section [V] the production and decay mechanisms of glueballs in hadronic 
reactions are discussed. In section PVTl the peculiarities of glueball production and behavior in Quark-Gluon Plasma 
(QGP) is considered. Section I VIII is dedicated to present two open topics the relation between the pomeron and 
glueballs and glueball-quarkonium mixing. Finally in section IVIIII we extract the main conclusions of our analysis 
and try to foresee future developments. 



II. LATTICE QCD 
A. Overview 



Gluon self-couplings in QCD suggest the existence of glueballs, bound states of mainly gluons. Investigating glueball 
physics requires an intimate knowledge of the confining QCD vacuum and it is well known that such properties 
cannot be obtained using standard perturbative techniques. To handle the nonperturbative regime of QCD, one 
can resort to numerical methods, known as lattice QCD. Lattice QCD needs as input the quark masses and an 
overall scale, conventionally given by Aqcd- Then any Green function can be evaluated by taking average of suitable 
combinations of lattice fields in vacuum samples. This allows masses and matrix elements, particularly those of weak 
or electromagnetic currents, to be studied. One limitation of the lattice approach is in exploring hadronic decays 
because the lattice, using Euclidean time, has no concept of asymptotic states. 

Lattice QCD was originally formulated by Wilson Q and is a clever implementation of the QCD dynamics using a 
finite difference formalism. The starting point is the correlation function in a discrete Euclidean space 

c(t) = <n|0+(t)0(o)|n> ~ [du [dip [ d^^2^( ^)^( x ^) e ' SFm ' SG(p ^ (!) 



where Sf is the fermion action and Sq the pure gauge action. The continuum limit is controlled by the input 
parameter (3 = 2N/g 2 (N is the number of color). By varying the inverse lattice coupling j3 we vary the lattice 
spacing a. The fermion action can be integrated out exactly in Eq. ([1]) to produce the fermion determinant. The 
determinant describes the dynamics of the sea quarks. In quenched QCD calculations, the determinant is set to a 
constant. (lr| 

The physics is extracted from the fit 

C(t) =^|(ft|0|n)| 2 exp(-M n i). (2) 

n 

\n) are the energy eigenstates and M n the corresponding masses. In practice, one has to choose the operator <fi which 
best overlaps with the lowest-lying glueball in the channel of interest. (f> is thus expanded in a basis with well-defined 
symmetry properties under the octahedral group and the variational coefficients are determined by Monte-Carlo 
simulations. [17| 

The spectrum in a box with periodic boundary conditions includes not only single glueball states, but also states 
consisting of several glueballs and torelons (gluon excitations which wrap around the toroidal lattice). Fortunately, 
torelons are found to overlap only weakly with single glueball states. [l7l Il8| 

Since (classical) gluodynamics is dimensionless, its observables will be also dimensionless. Masses are usually 
expressed in terms of the string tension m/y/a or the hadronic scale parameter r§m defined through the static 
potential between quarks [r 2 dV (r) / dr] r = ro = 1.65. [It], [H| Their values are very close and usually read 

v 7 ^ = 440 ± 20 MeV, r^ 1 = 410 ± 20 MeV. (3) 

In lattice calculations there are errors arising from the finite size of the lattice spacing a and the finite lattice volume. 
But for small enough a we expect the continuum limit to be approached as 

m(a) _ m(Q) , ^ 



with a constant c. 12011 
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FIG. 1: The mass spectrum of glueballs in pure SUc(3) gauge theory . The masses are given in units of the hadronic scale ro 
along the left vertical axis and in GeV along the right vertical axis. The mass uncertainties indicated by the vertical extents of 
the boxes do not include the uncertainty in ro. Numerical results are listed in the table below. In some cases, the spin-parity 
assignment for a state is not unique. The figure shows the smallest J value, the other possibilities are indicated in the second 
column of the table. 



B. Pure gauge spectrum 



The pur e gauge spectrum of quenched QCD was investigated initially by Morningstar and Peardon in an anisotropic 
lattice. [17| They used different spacings for the spatial a s and for the temporal at extensions with the ratio £ = a s /(H- 
This technique allows to control better the inherent errors induced by the lattice. 

Morningstar and Peardon identified 13 glueballs below 4 GeV. In order to distinguish single from multiple glueball 
states, they determined approximately the locations of the two-particles states using the mass estimation of the lowest 
few particles. In their estimates of these locations they assumed that the two glueballs do not interact and that the 
threshold for their production is given by the energy 

E 2G « \jp\ + m\ + yJpZ + ml, (5) 

where p\ = —p 2 and mi^ are the masses of the single glueballs. All states lying below the two-glueball threshold are 
then single glueball states. These authors pointed out that they cannot rule out a single glueball interpretation for 
higher states. They do not find any hint that their states are torelon pairs. 

Finite volume effects are quite under control. When going from a 6 3 x 40 lattice to a 8 3 x 40 lattice, the fractional 
changes on the mass 5 = 1 — ml /m are less than a few percent and consistent with zero. The largest effect of these 
errors occur in the and T-j*^ representations of the octahedral group. They are the main cause of uncertainties 

in the 0* ++ and 2* H glueballs. The proximity of the two glueball thresholds and the finite volume errors on the 
ylj ++ lead the authors to withhold judgment on whether or not this level is a single glueball. 

Lattice spacing errors, see Eq. are expected to be 0(a 2 , a*, a s a 2 s ) from perturbation theory. [l7l |2lj But the 
results for different £ suggest that 0(a 2 , a s a 2 s ) errors are negligible. They extrapolate to the continuum limit assuming 
0(a,g). In addition to these lattice errors of the dimensionless masses rom, one has to add the error arising from the 
scale parameters when presenting the absolute masses. 

Representations of the octahedral group are distinct from conventional spin representations of the Lorentz group. 
But one expects that in the continuum limit, the former match onto the latter. Once this extrapolation is achieved, one 
needs to identify the quantum number J PC of the lattice spectrum. The low-lying states do not lead to ambiguities. 
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TABLE I: Final continuum-limit glueball mass estimates mo- When a unique J interpretation for a state cannot be made, 
the other possibilities are indicated in the second column. States whose interpretation requires further study are indicated by 
a dagger. In column 3, the first error is the statistical uncertainty from the continuum-limit extrapolation and the second is 
the estimated uncertainty from the anisotropy. In the final column, the first error comes from the combined uncertainties in 
romc the second from the uncertainty in r^ 1 — 410(20) MeV. 





Other J 


TqTTIq 


m G (MeV) 


0++ 




4 21 (111(41 


1730 (50) (80) 


2++ 




5.85 (2)(6) 


2400 (25) (120) 


0"+ 




6.33 (7) (6) 


2590 (40) (130) 


0*++ 




6.50 (44)(7) f 


2670 (180) (130) 


1+- 




7.18 (4) (7) 


2940 (30) (140) 


2"+ 




7.55 (3) (8) 


3100 (30) (150) 






8.66 (4) (9) 


3550 (40) (170) 


o*-+ 




8.88 (H)(9) 


3640 (60) (180) 


3++ 


6,7,9,... 


8.99 (4) (9) 


3690 (40) (180) 


1 __ 


3,5,7,... 


9.40 (6) (9) 


3850 (50) (190) 


r -+ 


4,5,8,... 


9.50 (4)(9) f 


3890 (40) (190) 


2~~ 


3,5,7,... 


9.59 (4) (10) 


3930 (40) (190) 


3~ 


6,7,9,... 


10.06 (21) (10) 


4130 (90) (200) 


2+" 


5,7,11,... 


10.10 (7)(10) 


4140 (50) (200) 


o+- 


4,6,8,... 


11.57 (12) (12) 


4740 (70) (230) 



TABLE II: Glueball mass ratios. This ratios are not contaminated by anisotropy errors and are calculated using the empirical 
fact that correlations between different channels were found negligible. Note that the pseudoscallar glueball is clearly resolved 
to be heavier than the tensor. 

m(2 ++ )/m(0 ++ ) 1.39(4) 

m(0-+)/m(0 ++ ) 1.50(4) 
m(0* ++ )/m(0 ++ ) 1.54(11) 

m(l+-)/m(0 ++ ) 1.70(5) 

m(2-+)/m(0 ++ ) 1.79(5) 

m(3+-)/m(0 ++ ) 2.06(6) 
m(0*-+)/m(0 ++ ) 2.11(6) 
m(0~ + ) /m(2 ++ ) 1.081(12) 



The situation is different for higher states since they can belong to another multiplet with higher excitation not 
calculated. 

Their final results for the glueball spectrum are shown in Fig. [T] and in Tabic |TJ In the figure they assume the 
most likely spin interpretations. The table also contains any alternative spin attributions which cannot be ruled out. 
Several mass ratios are shown in Table HT1 which can be determined very accurately since they are not contaminated 
by anisotropy errors. Note that the pseudoscalar glueball is resolved to be heavier than the tensor. 

To convert lattice glueball masses into physical units the value of the hadronic scale tq must be specified. The 
estimate used, r^ 1 = 410(20) MeV, was obtained by combining Wilson action calculations of a/ro with values of the 
lattice spacing a determined using quenched simulation results of various physical quantities, such as masses of p and 
(j) mesons, the decay constant and the IP — IS* splittings in charmonium and bottomonium. A great deal of care 
should be taken in making direct comparisons with experiments since these values neglect the effects of light quarks 
and mixings with nearby conventional mesons. 

More recently Chen et al.[l9( have performed a similar calculation with larger lattices and larger volumes. We 
present in the Table [TTT1 a comparison of their results with those of Morningstar and Peardon.[l7j 

Meyer and Teper investigated also the pure gauge spectrum for (even) ++ states on a lattice in order to check 
the linearity of the Pomeron trajectory (see section |VII|) . [18] They reported also masses in other PC sectors. It is 
instructive to compare their results with the Morningstar and Peardon study. [17fl Although in these works absolute 
masses are expressed with different energy scales (rp 1 .fl7j and [3 ) i their close values [see Eq. (|3|)] allow the 
comparison of their absolute spectra displayed in Fig. [2] (left). Level orderings in both cases are identical but globally 
Meyer and Teper masses are smaller. Fig. [5] (right) presents mass ratios which are not contaminated by anisotropy 
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TABLE III: Continuum-limit glueball masses Ma for Chen et al. and for Morningstar and Peardon. 
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errors. In this case, errors bars are of the order of the symbols and are not shown. 
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FIG. 2: Comparison between Morningstar and Peardon (circles) and Meyer and Teper (triangles) of the mass spectrum of 
glueballs in pure SUc(3) gauge theory. Absolute masses (left) and mass ratios with respect the scalar glueballs (right). 

Recently, Meyer updated masses of the scalar and the tensor using his technique. [22| In this latter reference, he 
used the lattice scale ro allowing to compare with the other studies. The new masses are closer the Morningstar and 
Peardon's ones and read 

r Q M a ++ = 3.958(47), r Q M 2 ++ = 5.878(77). (6) 

All lattice calculations are now consistent and shown that, in pure gauge theory, the masses of the lowest states are 

M 0++ - 1.6 - 1.7 GeV, M 2 ++ - 2.4 GeV . (7) 

In full QCD interpolating operators for a state with given quantum numbers can also be constructed from quarks 
and anti-quarks. The pure glue operators might mix with the fermionic operators. If the mixing is very strong, 
the glueball masses obtained in this way, will have little to do with the glueball masses in the quenched calculation. 
Several methods have been applied to the interesting scalar sector, J PC = ++ , of the physical spectrum. 



Glueballs 



6 



Weingarten and Lee [23| studied the effect of the effect of quarkonium mixing with the glueball in the lowest ++ 
state in quenched QCD. The results were expressed as a mixing matrix 




where m g is the glueball mass, m a {s) is the mass of the ++ non-singlet ipip state, and E(s) is the mixing energy. 
Weingarten and Lee measured: m g = 1648(58) MeV, m a {s) = 1322(42) MeV, and E(s) = 61(58) MeV. The qualitative 
picture that emerges is that the /o (1710) is mostly ++ glueball, and the /o(1500) is mostly ss. The mixing energy 
E(s) has large lattice spacing errors. For example at a lattice spacing of a -1 ~ 1.2 GeV, the Weingarten and Lee 
result is E(s) ~ 0.36 GeV, while that of McNeile and Michael (24| is E(s) ~ 0.44 GeV. The analysis of Weingarten 
and Lee depends on the ++ states being well defined in quenched QCD. Bardeen et al.[25| have shown that there 
is a problem with the nonsinglet ++ correlator in quenched QCD. The problem can be understood using quenched 
chiral perturbation theory. The non-singlet ++ propagator contains an intermediate state of rj' — ir . The removal of 
fermion loops in quenched QCD has a big effect on the rj propagator. The result is that a ghost state contributes to 
the scalar correlator, that makes the expression in Eq. ([2]) inappropriate to extract masses from the calculation. 

A lattice QCD calculation that includes the dynamics of the sea quarks should reproduce the physical spectrum. 
Some insight into the composition of individual physical states, such as wether the physical particles couple to 
fermionic operators or pure operators could be studied as an effect of decreasing sea quark mass. Some studies have 
been performed for rif — 2 QCD [8j and it was found that the mass of the ++ glueball was reduced with respect 
to the quenched calculation by about 20%. Not so the tensor 2 ++ , whose value remained close to the quenched 
calculation. The mass of the ++ glueball on the UKQCD data set is degenerate with the mass of two pions.Q Due 
to the intricacies of the physical spectrum the lattice spacing used in the unquenched calculations must be reduced 
before direct contact can be made to phenomenology. 

In the real world glueballs have a decay width since they decay into two mesons. Lattice QCD calculations are 
performed in Euclidean space and this makes the computation of intrinsically complex quantities, as decay widths, 
complicated. (26| The decay width for the ++ glueball to two pseudoscalars has been calculated to be 108(28) MeV.[27j 
It is encouraging that it is small compared with its mass, however there is not yet consensus in this result. [28j 



III. CONSTITUENT MODELS 



Gluon self-coupling in QCD suggests the existence of glueballs. Incontrovertible experimental evidence for their 
existence remains elusive. A primary reason for this is the difficulty in extracting the properties of glueballs from 
the QCD lagrangian. We have seen that lattice QCD faces both computational and fundamental problems. We next 
describe a complementary way to describe glueballs, namely constituent models, which implement in a dynamical way 
the phenomenological properties of the confining QCD vacuum and the interaction among the gluons. 



A. The MIT bag model 



Hadrons are physical systems where quarks and gluons are confined in regions smaller than 1 fin. This experimental 
fact led the physicists of the MIT to develop a bag model of hadrons in QCD.[29j In their picture, quarks and gluons 
are confined in a bag, usually taken as a static spherical cavity. Confinement is described in the model by a boundary 
condition and a constant energy density B. The boundary condition makes the color flux vanish at the surface of the 
bag and it produces a quantization of the energy levels. B gives a constant energy term which stabilizes the bag at 
a finite size. Energy modes Ei = Xi/R are inversely proportional to the radius of the spherical cavity R. The energy 
in the cavity under these conditions, with rij massless constituents of type i, is 



AttBR 3 x - Xi 



R 

The bag energy |9j) encodes the masses of the states M but also the dynamics of the center-of-mass motion. One 
solution to this problem is the following: We consider that the bag is an eigenstate of the Hamiltonian H 2 = P 2 + M 2 . 
Then we simply take out the quantity 
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The mass then reads 

M 2 = E 2 - (P 2 ) . (11) 

In their application to hadrons,[30j the authors considered the bag constant as a free parameters. The minimization 
of the mass equation leads to a relation between the radius and B. Then the mass equation can be written in terms 
of only B and by fitting hadronic masses one can determine this parameter. 

Jaffe and Johnson were the first to apply this model to glueballs. [f| They found the modes of the gluon field in the 
cavity corresponding to the solution of the gluon equations of motion subject to the boundary conditions 

n M Gr = 0, (12) 
where n M is the normal to the bag surface and G% v is the gluon field strength. The two lowest modes are: 

Transverse Electric J p =1+ x TE = 2.744, (13) 

Transverse Magnetic J p =1" s T m = 4.493. (14) 

From these one obtains immediatly the masses of the low-lying states: (TE) 2 , ++ , 2++, M = 960 MeV ; (TE)(TM), 
Q-+, 2~+, M = 1.3 GeV ; (TE) 3 , 0++, 1+", 3 +_ , M = 1.45 GeV. Some authors pointed out the fact that lowest three- 
gluon glueballs ++ , F , 3 H have the opposite parity to that of the potential model predictions. [3lJ The parity of 
the TE mode 1 + causes this difference. The lattice spectrum seems to support the bag model. 

Differently to potential models, constituents in the bag model are not confined by a potential. Particles are almost 
free inside the bag and the confinement is ensured by the boundary conditions (|12p . A more elaborated description 
should lift degeneracies. Carlson et al. [32| added to the bag energy ([9]) a correction AE representing the spin 
splitting induced by one-gluon-exchange interaction. This shift includes tree-level scattering diagrams but also self- 
energy contributions, 

AE = Y, AE ij + ^2 AE h ( 15a ) 
AEij = -^-^(ti-tj(a,ijSi- Sj + bijTij+CijIij) + dijVij), (15b) 

H 

AE, = -*&-t 2 e h (15c) 
H 

Here tj and Si are the generators of color and spin respectively. V is the projector onto the color-octet spin-one state. 
J is the identity operator and a(R) is the running coupling constant, 



. „. 27T7J 1 

v ' 9 ln[l + l/(Ai?)«] v ' 

This Ansatz for a(R) simulates the saturation for large R. n is a positive integer parameter. The authors used n = 2 
and A = 0.172 GeV in their calculations. The coefficients a, b, c, d, e in the relations (|15[) are given in the reference. 

In the original version of the bag model, the bag constant was a free parameter determined from data. Hansson 
et al. proposed a model for the QCD vacuum wich allows B to be calculated given a(R) and exE (the self-energy of 
the lowest TE gluon mode).[33j The basic idea is that the vacuum is filled with ++ (TE) 2 glueballs which form a 
negative energy condensate. The expression for B 

B =s4 ( - mJ)1/a (17) 

relates its value to the (negative) mass squared of the scalar glueballs m 2 . The excitation of this condensate gives rise 
to a observable scalar glueballs. [34{ Their results are shown in Fig. [3j This figure displays also the glueball spectrum 
obtained in the bag model of Chanowitz and Sharpe along similar lines. [35[ Masses in the bag model are globally 
lower than in lattice QCD. 

In bag models, particles are confined in a spherical bag. This approximation is only valid for totally symmetric 
,7 = states. Moreover, spherical glueballs are not stable since a vector field can never give rise to a spherically 
symmetric pressure. Robson pointed out this flaw and developed a toroidal bag model for glueballs. [3^] This picture 
of glueballs is close to the flux tube model of Isgur and Paton. [37| 
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FIG. 3: The mass spectrum of glueballs in bag models: Jaffe and Johnson 6 (triangles), Carlson et al. 32 (blue circles) and 
Chanowitz and Sharpe 35 (squares). 



B. Gluon mass 



A possible way to handle glueballs is to consider massive quasi-gluons interacting via a QCD inspired dynamics. 
The gauge bosons are massless at the Lagrangian level but there are hints that they obey massive dispersion relations. 
The gluons remains massless to all orders in perturbation theory. However, nonperturbative effects, like confinement, 
and their self-interactions, can be described by a constituent mass. Several definitions exist for this constituent mass. 

The so called dynamical mass, is defined by the position of the pole of the dressed gluon propagator. Cornwall 
arrived to such a dynamical mass by analyzing the gluon Dyson-Schwinger equations in the early 80's.[38| This 
infinite set of couple equations cannot be solved analytically. One must resort to a truncation scheme. By a clever 
resummation of Feynman diagrams, Cornwall found a gauge-invariant procedure to deal with these equations. With 
this technique, a full gluon propagator in quarkless QCD emerges 



d-V) = (<7 2 + m 2 (gWln 



with a dynamical mass 



2/ 2\ 2 

m (q ) = m 



'ln[(g 2 +4m 2 )/A 2 
In (4m 2 /A 2 ) 



■ 4m 2 (q 2 ) 
~ 2 



-12/11 



(18) 



(19) 



In Eq. (p~8|> . b — 11A/487T 2 is the first coefficient of the beta function for quarkless QCD and the mass term that 
appears has finite value at zero momentum. The gluon mass can be related to the gluon condensate (G* G£") from 
which the value m = (500 ± 200) MeV arises. 

Bernard proposed a different definition for the gluon mass.[39| Consider the potential energy of a pair of heavy, 
static sources in the adjoint representation of the color group. As the separation of the adjoint sources (static gluons) 
is increased, the potential will increase linearly as a string or a flux tube is formed between them. The energy stored 
in the string will at some point be large enough to pop up a pair of dynamical gluons out of the vacuum. The effective 
gluon mass is defined as half of the energy stored in the flux tube at this point. Monte-Carlo simulations of this 
phenomenon show a effective gluon mass in the range 500-800 MeV. 

The effective gluon mass was also investigated in the bag model. [4f| Even though the gluon is massless in the bag 
model, a net energy of 740 ± 100 MeV is required to produce a gluon due to confinement. 

All these arguments support the use of an effective gluon mass to describe the dynamics of QCD. It is therefore 
possible to envisage an approach to bound states made of constituent massive gluons. Since two-gluon glueballs have 
always a positive conjugation charge, a study of the full spectrum must include also three-gluon glueballs. 
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C. Two-gluon glueballs 



One of the pioneering works on two-gluon glueballs was the study by Cornwall and Soni.,4l| The large value of the 
effective gluon mass led them to propose a nonrelativistic approach to gluonium. They used a confining potential 
which saturates at large distances constrained by Bernard's results 39] 

V c (r) = 2m(l-e- r / r ^ . (20) 

This screened potential led to a spectrum with relatively low glueball masses. The scalar and tensor glueballs had 
masses around ~ 1.3 and ~ 1.6 GeV, respectively. They used for the string tension twice the value commonly 
used for mesons, a decision which they justified by arguing that a gluon acts as a qq pair. The Coulomb and 
spin-dependent interactions at short-range were derived from a nonrelativistic expansion of the Feynman graphs for 
two-gluon scattering. They considered massive exchanged gluons (with the same mass as the constituent one) to keep 
the gauge invariance of the amplitudes to the given order. This one-gluon exchange (OGE) potential involves Yukawa, 
spin-orbit, spin-spin and tensor forces, 



Kge(r) 



2s - 7m 2 



3A 
2^2 



L-S- 



6m 2 
e~ 



I* 

A 



r dr 



2m 2 



Tr\6(r) /4m 2 -2s 
3m 2 V m 2 



(S- V) 2 -^5 2 V 2 



(21) 



s is the glueball mass squared, which we can set to s = 4m 2 in a first approximation, and A = 3a s is the strong 
coupling constant. 

Cornwall and Soni presented results for states with quantum numbers L — 0, J PC = ++ ,2 ++ , and L = 1, 
J PC = 1 ^,2 Despite the fact that the gluons acquire a mass they remain transverse. For transverse particles 
the J = 1 states are forbidden as a consequence of Yang's theorem, [42[ i.e., a gluon in a massless representation has 
only two projections for its spin and therefore two transverse gluons cannot bind into a J = 1 state. Thus we must 
incorporate this feature into the above formalism. 

When dealing with massless representations, the conventional J = L + S decomposition is not useful anymore. 
The formalism to treat two-body relativistic scattering developed by Jacob and Wick[43] allows also the description 
of representations with only transverse gluons. We sketch its main features and apply this formalism to the study of 
the two-gluon glueball. 

The Jacob and Wick formalism is based on states, |J,M;Ai,A 2 ), which are eigenstates of J 2 and J z and where 
Ai and A2 represent the two allowed spin projections. In the case under consideration, the projections can only be 
maximal, i.e. ±s for a particle with helicity-s . The angular part of these states are related with the conventional 
basis states by means of Clebsch-Gordan coefficients: 



|J,M;A l! A 2 )=^ 



L.S 



2L + 1 



-1 1/2 



2J+ 1 



{LS0A\ J A ) («i« a Ai(-Aa)| SA ) \ 2S+1 L.,) , 



(22) 



with A = Ai — A2. The radial part (depending on J) is determined variationally with the Hamiltonian. The states 
are not eigenstates of parity. For a two-gluon state, s\ = S2 = 1, it holds 



P|J,M;A 1 ,A 2 ) = (-l) J |J,Af;-A 1 ,-A 2 ). (23) 

The key point is that physical states must not only be eigenstates of the total angular momentum operator but they 
must also be eigenstates of parity. Such a requirement is fulfilled by the following linear combinations 

\H±; J p ; A x , A 2 ) = i= { \J, M; A x , A 2 > ± \ J, M; -A x , -A 2 )} , (24) 

for which P \H±; J p ; X u A 2 ) = P \H±; J p ; Ai, A 2 ), with P = ±(-l)' 7 . In the latter, the \H±; J p ; Ai, A 2 ) states will 
be referred as helicity states. When the two particles are identical, the wave function should be an eigenvector of the 
permutation operator Pi 2 . The basic state (|24p are eigenstate of the permutation in the case of massless particles 
since Ai = ±A2 and[43[ 



P 12 |J,M;A 1 ,A 2 > = {-1) J \J,M;X 27 \ 1 ) 



(25) 
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A system of two gluons has to be totaly symmetric and this constraint leads to selection rules on the spin and the 
parity. Indeed, for singlet states (defined by Ai = A 2 = A), we have 

Pia \H±; J p ; A, A) = (-1) J \H±; J p ; A, A) , (26) 

leading to the families with the even spin and positive or negative parity . The minimum spin in this case is J = 0. 
For the doublet (defined by Ai = — A 2 = A), the minimum spin is J — 2 and we have 

Pia \H ± ; J p ; A, -A) = ±(-l) J |if±; J p ; A, -A) . (27) 

This time we can have odd spin (the lowest is 3 + ) but also even spin with positive parity. We thus observe the 
emergence of selection rules according to the value of J and P. 

The results for states of two transverse gluons are formally identical to the case of states made of two photons, i.e., 
Yang's theorem. [42j The total color wave function is assumed to be a singlet, which is totally symmetric, and does not 
explicitly appear in the computations. Taking into account the Bose symmetry, one finds that there are four allowed 
helicity states, namely 

\s ± ;j p ) = \H±-,j p ) Xa=Xi , ( 28a ) 

\D ± ;J p ) = \H ± :J p ) X2 _ Xi . (28b) 

The selection rules impose restrictions to the possible values of the total angular momentum and parity of these four 
types of states. It can be checked that one can only obtain the following states 

\S+;(2k)+), |S_;(2fc)-), \D+; (2k + 2)+) , |D_; (2k + 3)+) , k e N. (29) 

The S- and D-labels stand for helicity-singlet and -doublet respectively. We recognized in Eq. (f2"9")) the four families 
predicted by Yang. [52] 

It should be noted that a state made of two gluons in a color singlet state has always positive charge conjugation 
(C = +1). More explicitly, the states in Eq. (|29|) give rise to the following glueball states 

\S+;(2k)+) => ++ ,2 ++ ,4 ++ ,... (30a) 

\S- ; (2ky ) => (T + ,2-+4- + ,... (30b) 

\D + ;(2k + 2) + ) => 2++,4 ++ ,... (30c) 

|L>_;(2fc + 3) + ) => 3 ++ ,5++,... (30d) 

It is readily observed that only the \S±; (2k) + ) states can lead to J = 0, while the \D±) states always have J > 2 
(since J > |Ai — A2I). Obviously, a consequence of Yang's theorem is that no J = 1 states are present. Only the \D-) 
states can generate an odd- J, but J is at least 3 in this case. 

Lattice QCD confirms the absence of the 1 ^ and 1 ++ states, at least below 4 GeV. It is worth mentioning that 
glueball states with even-J and positive parity can be built either from the helicity-singlet or from the helicity- 
doublet states. The important result is that the gluons remain transverse and therefore the helicity formalism exactly 
reproduces the J PC content for glueballs which is observed in lattice QCD, without the extra states which are usually 
present in potential models. 

The helicity formalism was applied for the first time by Barnes. [44| It has several advantages not shared by the more 
conventional nonrelativistic LS'-basis. It avoids spurious states forbidden by the coupling of two transverse gluons 
but also reproduces the lattice QCD hierarchy, i.e. 

±+ , 2++, 2 ±+ , 3++, 4++, 4 ±+ , 5++, 6++. (31) 



Within this approach, a given J PC state can be expressed as a linear combination of (L, S) states thanks to Eq. 
The complete expressions for these decompositions can be found in Mathieu et aL[45( We give here the angular 
dependence of the ground states of Eq. (|3H)) : 



|S+;(0)+) = V^| l5 o) + V^| 5 ^o), (32a) 

\S-;(0)-) - | 3 P ), (32b) 

\D + ;(2) + ) = y||5s 2 ) + y||^D 2 ) + yi| 5 G 2 ), (32c) 

|^-;(3) + ) = yfr^3> + yf | 5 G 3 ). (32d) 
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These decompositions are essential for computing the matrix elements of relativistic operators (spin-spin, spin-orbit 
and tensor). Let us note that the matrix elements of these operators are equal for \S + ; (2k) + ) and (2fc)~).[45| 

Even though in this approach the singlet states J p = (2fc + ,2fc~) are degenerate, with a Cornell- type (linear + 
Coulomb) potential, a nonrelativistic kinetic energy, which incorporates an ad hoc gluon mass to, and using the helicity 
formalism, Barnes was able to reproduce the qualitative feature of the pure gauge sector finding M(0 ±+ ) = 4.36 to. 
The higher mass ratios were not in perfect agreement with modern lattice results, implying the need for modifications. 

We emphasize that considering transverse gluons is essential for finding the correct hierarchy. However, transverse 
particles are generally massless and even if nonperturbative effects are able to give a mass to the gluon, one may 
wonder if a nonrelativistic kinetic energy p 2 /2m is appropriate. Indeed, the nonrelativistic kinetic energy is just the 
limit of the more general Dirac operator \/p 2 + to 2 . This semi-relativistic energy is also valid for massless particles 
such as gluons. 

Brau and Semay compared different models for glueballs. [46[ Models with a nonrelativistic kinetic energy were not 
able to reproduce correctly the lattice gauge spectrum for realistic values of the parameters. They concluded that 
a semi-relativistic Hamiltonian, i.e. 2y / p 2 _|_ V(r), is an essential ingredient to handle glue states. Nevertheless, 
all models analyzed used an LS'-basis and were plagued with unwanted states. Moreover, they had to include the 
short-range potential Eq. (|21[) to lift some degeneracies between the states. This is not the case if we implement 
transverse gluons by means of the Jacob and Wick formalism, then automatically the degeneracies are lifted. 

This improvement was carried out in a work based on the Coulomb gauge Hamiltonian where a relativistic kinetic 
energy was used.[47j In this model, gluons are linked by an adjoint string. The adjoint string tension a a = (9/4) cr is 
expressed in terms of the well-known fundamental string tension for mesons a through the Casimir scaling hypothesis 
supported by lattice calculations. [HI Using typical values for the parameters, a = 0.18 GeV 2 for the fundamental 
string tension (extracted from mesons Regge trajectory) and as = 0.4 for the strong coupling, this model encodes 
the essential features of glueballs. 

The spectrum of the Coulomb gauge Hamiltonian was in good agreement with lattice QCD. Moreover, the singlet 
2 h and 2 ++ are degenerate as in the Barnes' model, a characteristic of the helicity formalism. The authors found 
a difference between the scalar and pseudoscalar glueballs. This splitting, about 250 MeV was nevertheless not as 
strong as in lattice QCD (800 MeV). 

Recently, this problem was revisited keeping the basic ingredients needed for obtaining an acceptable pure gauge 
spectrum compatible with lattice results, i.e. semi-relativistic energy and the helicity formalism for two transverse 
gluons. [45[ A simple Cornell potential was used but an instanton induced force was added and with it the splitting 
between the scalar and pseudoscalar glueballs was reproduced. There are arguments favoring an attractive (repulsive) 
interaction induced by instantons in the scalar (pseudoscalar) channel of glueballs. [ii^] 

In all the constituent models for glueballs, the confining interaction follows from phenomenological considerations 
(breakable strings or linear potentials) . But it is also possible to derive an effective Hamiltonian for bound states from 
the QCD Lagrangian. Kaidalov and Simonov used the field correlator method to extract a relativistic Hamiltonian 
which describes an adjoint string with gluons at its ends. [5(| The introduction of auxiliary fields (or einbein field) /j, 
and v leads to 



L(L + 1) 



1)2 
2> 



vdf3 



1 a 2 A d[3 



2;/ 



ud/3. 



(33) 



In Eq. (|33| . /i is introduced to get rid of the square root and (fi) = {yfp*) represents the effective gluon energy, (y) is 
the energy stored in the rotating string. The einbeins are eliminated using the equations of motion, S^Hq = S u Hq = 0. 
For L = 0, we find v = a at, and the replacement of the auxiliary field leads to 



H Q = 2vV + a A r. (34) 

For L ^ it is not possible to eliminate v analytically and the optimization S^Hq = should be done numerically. 

In this approach, gluons are massless but gain an effective mass fj,o = (fi) ~ 0.5-1 GeV. The expectation value 
is taken on an eigenstate of Hq therefore fio is state-dependent. The spin splitting operators are corrections of 
order ^ 2 and are computed perturbatively. These corrections contain the conventional structure arising from the 
one-gluon-exchange and, in addition, a spin-orbit term coming from the Thomas precession, 
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AH 



(35) 



A A L S 87rA r . , a a 

+ -^(3(Si-f)(S 2 -r)-S 1 -S 2 ) 
r Mo 

with A = 3a s . 

The numerical results for the two-gluon glueballs in all the models discussed are displayed in Table HVl and Fig. 2] 
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FIG. 4: Comparison between lattice results 17 (circles) and two-gluon glueballs (left) from Ref. 47 (squares) and Ref. 45 (tri- 
angles) and three-gluon glueballs (right) from Ref. 55 (squares) and Ref. 57 (triangles). 



TABLE IV: Two-gluon glueballs spectra in different models. In brackets, the ratios with respect to the scalar glueball. Masses 
are in GeV. 





ref. [44] 


ref. [471 


ref. [45] 


ref. [511 


ref. [17, 19] 


ref. [18] 


0++ 


(1.00) 


1.98(1.00) 


1.72(1.00) 


1.41(1.00) 


1.71(1.00) 


1.48(1.00) 




(1.40) 


3.26(1.65) 


2.54(1.48) 


2.41(1.71) 


2.67(1.56) 


2.76(1.89) 


o-+ 


(1.00) 


2.22(1.12) 


2.62(1.52) 


2.28(1.62) 


2.56(1.50) 


2.25(1.54) 




(1.40) 


3.43(1.73) 


3.44(2.00) 


3.35(2.38) 


3.64(2.13) 


3.37(2.31) 


2++ 


(1.16) 


2.42(1.22) 


2.59(1.51) 


2.30(1.63) 


2.39(1.40) 


2.15(1.47) 




(1.35) 


3.11(1.57) 


3.08(1.79) 


3.32(2.35) 




2.88(1.97) 


2" + 


(1.35) 


3.09(1.56) 


3.08(1.79) 


2.70(1.91) 


3.04(1.78) 


2.78(1.90) 




(1.68) 


4.13(2.09) 


3.73(2.17) 


3.73(2.65) 


3.89(2.27) 


3.48(2.38) 


3++ 


(1.42) 


3.33(1.68) 


3.25(1.89) 




3.67(2.15) 


3.39(2.32) 


4 ++ 


(1.63) 


3.99(2.02) 


3.77(2.19) 






3.64(2.49) 




(1.71) 


4.28(2.16) 


3.96(2.30) 








4"+ 


(1.71) 


4.27(2.16) 


3.96(2.30) 








5++ 






4.21(2.45) 








6++ 






4.60(2.67) 






4.36(2.98) 



The low-lying positive C-parity glueball states, seem to favor a constituent picture with two gluons interacting via 
a linear potential, i.e. linked by a string. A relativistic approach with two transverse gluons leads to a spectrum in 
good agreement with the lattice hierarchy. 

This stringy picture leads to a Regge trajectory 



2-KG 



(36) 
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well-known in the meson sector. The experimental slope for mesons \/(2ira) « 0.9 GeV~ 2 corresponds to the typical 
value for a fundamental string tension a ss 0.18 GeV 2 . When we extend this picture to glueballs, the slope rises to 
0.4 GeV~ 2 if the Casimir scaling hypothesis is used. If one argues that even spin positive C-parity glueballs lie on 
the pomeron trajectory, a problem arises. Indeed, the experimental soft pomeron slope is a 1 ~ 0.25 GeV -2 . Since the 
pomeron trajectory carries physical particles, a first solution is that glueballs probably mix with quark states. Let us 
note before finishing this discussion that to describe negative C-parity states in constituent models one needs at least 
three gluons. 

A very different approach to the interpretation of the spectrum is in terms of a closed flux tubes. This picture was 
introduced by Isgur and Paton. [37| This model is composed by a closed loop of fundamental flux with no constituents 
gluons at all. Such a loop has phononic and orbital degrees of excitation. These two modes lead to two Regge 
behaviors at large spin, 



phononic : J ~- M , 

8lT(T 

In either case one obtains a slope 0.2-0.3 GeV~ 2 which is in the right range for the pomeron [see Eq. (162p ].(l8l] 
Another interesting feature of the closed flux tube model is the appearance of low-lying odd spin PC = H — states in 
the spectrum. In their original paper. Isgur and Paton discussed masses for some low-lying pure glue states (see Fig. 
[Vp .[37j|. These authors argued that the true values of their parameters were not known and they simply chose them 
to fit expected results. The splittings are not sensitive to their choice of parameters but their absolute values are. 



orbital : J 



3V3 



M 2 . 



(37) 



TABLE V: Low lying glueball states in the flux tube model. 



J™ Mass (GeV) 

IF* L52 

l +_ 2.25 

0++ 2.75 

++ ,0+-,0- + ,0— 2.79 

2 ++ 2.84 

2 ++ ,2 ++ ,2 ++ ,2 ++ 2.84 

1+- 3.25 

3+" 3.35 



D. Three-gluon glueballs 

A complete investigation of the glueball spectrum in constituent models has to include three-gluon glueballs. Indeed, 
in this approach negative C-parity glueballs involve at least three constituents. There are two color wave functions, 
totally symmetric or totally antisymmetric, coupling three adjoint representations into a singlet. They do not mix 
and we are only interested in the symmetric one, d a bcA a ^A b v A c p , namely the C = — states. 

Soni with Hou extended his paper with Cornwall to three massive gluons. [HJ Their potential is a generalization of 
the previous one, 



V 999 = E + ( 38 ) 

i<j 

with Vq the confining potential of the breakable string of Eq. (|2TJ)) . A factor one-half is added because one needs to 
remove three (and not six) gluons from the vacuum to screen the three gluons that are originally there in the glueball. 
The short-range potential is not the OGE potential Eq. ([2Tj) with a different color factor, because one should also take 
into account the annihilation diagram not present in two-gluon glueballs, but it involves the same structures (spin- 
spin, tensor and spin-orbit interactions). The nonrelativistic kinetic energy and the screened confining potential lead 
to a spectrum which is too low. The low-lying states are nearly degenerate with a mass 4.8 times the constituent mass 
m of the gluons. By low-lying, the authors mean that every pair of particles is an S'-wave. Symmetry considerations 
on the total wave function imply that the low- lying states are the h , 1~~, 3~~. 
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The 3 is the lowest maximum spin state. One often argues that such states with maximal spin lie on the odderon 
trajectory, the counterpart to the pomeron trajectory in hadron-hadron scattering at high energy. The linearity of 
the odderon trajectory was checked using a Coulomb gauge Hamiltonian. [55j] Their masses for the 1 , 3 , 5 , 
are shown in Table IVIl 



TABLE VI: Odderon quantum numbers and masses in MeV. 





Conf. 


a (GeV 2 ) 


1— 3"" 


5— 


7— 


s 






1 3 


3 


3 


L 









2 


4 


Llanes-Estrada[55] 


A 


0.18 


3950 4150 


5050 


5900 


Simonov [50J 


A 


0.238 


3490 4030 






Simonov[511 


A 


0.18 


3020 3490 


4180 4960 


Y 


0.18 


3320 3830 


4590 


5250 


Mathieu[56, 571 


Y 


0.21 


3999 4167 


5263 




Latticefl9l 




0.1939 


3830 4200 






Lattice[48l 




0.1939 


3100 4150 







In their Coulomb gauge approach to three-gluon glueballs, these authors chose a spin-independent potential which 
is a sum of two-body Cornell ones and a A-shape for confinement. Another Ansatz for the confinement, the K-shape, 
is sometimes preferred. The y-shape is the generalization of confinement in baryons where every quark provides 
a flux tube. These flux tubes coming from gluons meet in a point where the total energy (or the length since the 
energy density is constant) is minimal. It is worth mentioning that, under the Casimir scaling hypothesis, a simple 
application of the triangular inequalities shows that the A-shape is energetically more favorable than the y-shape.[52| 
This demonstration was recently confirmed by a lattice study of the three-gluon potential. [53| 

Kaidalov and Simonov investigated both Ansatze for confinement with a nonrelativistic Hamiltonian (5l| 

H ggg = P * + ri +P * + f + V A , Y (ri,r 2 ,r 3 ), (39) 

with 

9 3 

V A (ri,r 2 ,r 3 ) = ct^ \r t - r 3 \, V Y (r 1 ,r 2 ,r 3 ) = -a^2\r t - R Y \. (40) 



4 

Kj 1=1 



They found the eigenvalues of this operators by using the hyperspherical formalism. The spin-average masses of this 
Hamiltonian for the odderon states (odd) are displayed in Table PVTl Clearly, the y-shape potential leads to a 
higher spectrum than the A-shape one. 



TABLE VII: Glueball spectrum in the PC = H — sector. 





Latticcfl9] Mathicu[57] 


Simonovf51] 


o+- 


4780 4656 


4090 


1+- 


2980 4626 


4090 


2+" 


4230 4542 


4090 


3+" 


3600 4568 


4090 


5+- 


4110 [18] 5317 





As we saw in the case of two-gluon glueballs, a nonrelativistic description of such states is not appropriate. However, 
in the negative C-parity sector, the lack of lattice results for high spin states does no allow us to draw any definitive 
conclusions. Also, the Regge trajectory of the odderon is not well understood yet. New theoretical and experimental 
research on these topics would be helpful for understanding these glueballs. 

On the other hand, lattice QCD exhibits an interesting spectrum in the PC = H — sector. There is nearly 2 GeV 
between the lowest l" 1 state and the highest H . According to the LS decomposition in constituent models, all 
these states (0 -1 , 2 H ,3 H ) should be L = 1 and degenerate. This feature appears clearly in Simonov and 
Kaidalov's paper. [511] 
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Using a generalization of the Hamiltonian Eq. (|39|) but with a semi-relativistic kinetic energy, Mathieu et. al found 
also a degeneracy between the J H states in disagreement with the lattice results. [57] Their Hamiltonian involves a 
Y-junction at the center of mass for the confinement and an OGE potential for the short-range part 



3 3 

H = Y,ytf + a Y,\ r *- R °\ + v ° G v- (41) 

i=l i=l 

Their study is an extension of a first one where only L — states were considered. [56jj The parameters for this 
three-gluon glueball model were fitted on the two-gluon glueballs of a previous work.[46j Previously, the 2~~ had been 
computed with the wrong wave function. 15511 In this study, the right 2 state was obtained. It was noted that this 
result disagrees with that of lattice QCD.[56[ In a model with longitudinal gluons, the 2 cannot lie between the 1~~ 
and the 3 __ as lattice calculations have shown (cf. Figd]). This was a hint that a description with the nonrelativistic 
decomposition J = L + S should be inadequate to handle the pure gauge spectrum. In the PC = H — sector, all 
states within this approach lie in the same range and they contradict the lattice results. Hence, we conclude that a 
description of many-gluon states with a LS'-basis is not appropriate. A solution could be to implement the hclicity 
formalism for three particle states for three transverse gluons. One hopes that then it would be possible to reproduce 
the correct hierarchy. Indeed, the lowest states with three longitudinal gluons are the 1 and 3 while in lattice 
QCD the lowest negative C-parity are the 1 H and 3 H . 



E. AdS/QCD 

An alternative approach to strong interactions is based on the idea that they have a description in terms of 
strings. A remarkable step in this direction was given by Maldacena proposing the equivalence between conformal 
fields and string theory in anti-de Sitter spacetime (AdS/CFT correspondence) . [58| In particular glueball operators 
of the conformal gauge theory defined on the AdS boundary are in correspondence with the string dilaton field. The 
description of the strong interactions based on this correspondence requires the breaking of conformal invariance, 
which can be done in different ways. The spectrum of the glueball has been calculated in two of these models for 
breaking. The first introduces a Schwarzschild black hole in AdS to break scale invariance. [59| The corresponding 
supergravity equations do not admit closed form solutions and the spectrum has been calculated using approximate 
methods (see Table IVIIip . [gO, HH The second possibility is to imitate the Randall-Sundrum model by considering two 
slices of AdS sticked together assuming that there is a bulk/boundary correspondence between the glueballs and the 
dilaton. [62j Imposing boundary conditions the string dilaton acquires discrete modes and these modes become the 
spectrum of glueballs (see Table IVIII[) . 

Let us for the purposes of illustration, describe the analysis of the complete glue ball spectrum calculation for the 



AdS 7 black hole supergravity dual of QCD4 in strong coupling limit: g 2 N — > 00. 61] Despite the expected limitations 
of a leading order strong coupling approximation, the pattern of spins, parities and mass inequalities bare a striking 
resemblance to the known QCD4 glueball spectrum as determined by lattice simulations at weak coupling. 



To approach QCD4 one begins with M theory on AdS 7 x S 4 . One compactifies the "eleventh" dimension (on a 
circle of radius R\ ) to reduce the theory to type IIA string theory and then following the suggestion of Witten raise 
the "temperature", /3 , with a second compact radius R2 in a direction r, with j3 — 27ri?2- On the second "thermal" 
circle, the fermionic modes have anti-periodic boundary conditions breaking conformal and all SUSY symmetries. 
This lifts the fermionic masses and also the scalar masses, through quantum corrections. The 't Hooft coupling is 
g 2 N — 2Trg s Nl s / R2, in terms of the closed string coupling, g s and the string length, l s . Therefore, in the scaling 
limit, g 2 N — > 0, if all goes as conjectured, there should be a fixed point mapping type IIA string theory onto SU(N) 
pure Yang-Mills theory. 

One considers the strong coupling limit at large N, where the string theory becomes classical gravity in the AdS 
black hole metric, 

ds 2 = (r 2 - ±)dr 2 + r 2 ]T dx 2 l+ (r 2 - l^dr 2 + ±dQ 2 , (42) 

i=l,2,3,4,ll 

with radius of curvature, R\ds = 87n? s A7;?. The dimensionful parameters have been removed in the metric by a 
normalization setting Ra<is = 1 and P — 27r/3. 
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TABLE VIII: Masses of glueball states J PC with even J expressed in GeV, estimated using the sliced AdSs x S 5 space with 
Dirichlet (Neumann) boundary conditions. 63 The mass of ++ is an input from lattice results. 17 ' 64 



Dirichlet 


lightest 


1 st excited 


2 nd excited 


glueballs 


state 


state 


state 


0++ 


1.63 


2.67 


3.69 


2++ 


2.41 


3.51 


4.56 


4++ 


3.15 


4.31 


5.40 


6++ 


3.88 


5.85 


6.21 


8++ 


4.59 


5.85 


7.00 


10++ 


5.30 


6.60 


7.77 



Neumann 


lightest 


1 st excited 


2 nd excited 


glueballs 


state 


state 


state 


0++ 


1.63 


2.98 


4.33 


2++ 


2.54 


4.06 


5.47 


4 ++ 


3.45 
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FIG. 5: The AdS glueball spectrum 61 for QCD4 in strong coupling (left) compared with the lattice spectrum 17 for pure SU(3) 
QCD (right). The AdS cut-off scale is adjusted to set the lowest 2 ++ tensor state to the lattice results in units of the hadronic 
scale 1/ro = 410 Mev. 



The strong coupling glueball calculation consists of finding the normal modes for the bosonic components of the 
supergraviton multiplet in the AdS 7 x S 4 black hole background. We are only interested in excitations that lie in the 
superselection sector for QCD±. Imposing this restriction and exploiting symmetries of the background metric reduce 
the problem to six independent wave equations which have to be solved with the appropriate boundary conditions. 

There is a rather remarkable correspondence of the overall mass and spin structure between strong coupling glueball 
spectrum and the lattice results at weak coupling for QCD4, see Fig. [5] Apparently the spin structure of type IIA 
supergravity does resemble the low mass glueball spin splitting. The correspondence is sufficient to suggest that the 
Maldacena duality conjecture may well be correct and that further efforts to go beyond strong coupling are worthy 
of sustained effort. 

The basic idea behind the AdS/CFT correspondence is that the low mass glueball spectrum can be qualitatively 
understood in terms of local gluon interpolating operators of minimal dimension an idea which is not new. [11, l65j 
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These operators are in rough correspondence with all the low mass glueballs states, as computed in a constituent 
gluon or bag models discussed previously. Thus the AdS/QCD if valid would justify the behavior of the gluon as a 
constituent. Let us note, however, that in constituent models, confinement is ensured by the potential and in the bag 
model and AdS/QCD, by a boundary conditions. 



IV. GLUEBALLS AND QCD SUM RULES 



The QCD sum rule (SR) approach is one of the most widely used methods to obtain the information about glueball 
properties. (6^, [6?], [fil, [69| It is based on the operator product expansion (OPE) of the correlator of two glueball 
interpolating currents J G 

n(Q 2 ) = * / ' d A x e lq - x (0\TJ G (x)J G (0)\0) (43) 

in the deeply-Euclidean domain Q 2 = —q 2 ^> A.q GD . The method takes into account perturbative, as well as, 
nonperturbative gluonic contributions to this correlator. The perturbative contributions arise by direct calculation 
of Feynman diagrams, while the non perturbative contributions are associated, as we shall discuss, with vacuum 
expectations values of the correlators, i.e. condensates, and sometimes with direct instanton contributions. 

The interpolating currents, consistent with the minimum number of gluon fields, used to study low mass glueballs 
are the field stre ngth squared (S = ++ ), the topological charge density (0 h ) and the energy- momentum tensors 
(T = 2++): [H3 

Js(x) = a s G^(x)G^(x), (44) 
Jp(x) = a s G° v (x)G* v (x), (45) 



JP{x) = Gl a {x)G a va {x) + 9 —G% a {x)G% a {x), (46) 



4 



where G^x) = (l/2)e^ af3 G a afj (x). 
The dispersion relation 



^_ (-Q 2 r r° . imn(s) 

l %"( S + Q 2 ) 



n(Q 2 )-^^/ ds _ n ,_ y> +a + fll Q 2 + ..., (47) 



where a.; are substraction constants, allows us to connect the theoretical side of the SR (left hand side) with observed 
properties of the glueballs introduced through ImII(s) (right hand side). Conventionally, for the phenomenological 
part of the SR, the following form for the spectral representation is used, 



Imn(s) phcn = Kf Gi m G .5(s - m G z ) + - s o) Imn(s) thoor , (48) 

i 

which corresponds to a sum over narrow width resonances located at Sj = M Gi plus a continuum at large energy 
s > So- In Eq. [j4"8"]). fa 4 is the residue of the i th -glueball state defined by the following matrix element, 

(0\J G (0)\G i )=m%J Gi . (49) 

Here, m Gi is i th -glueball mass, so is the continuum threshold and II(s) thcor is the perturbative part of the correlator. 
This is, of course, a rather simplified model for the physical spectral density. However, it has been shown that such 
a model gives a rather good description of the mass spectrum for the ordinary hadronic states. One may expect 
therefore, that such a model can be also used for the extracting of the glueball masses. 

Before carrying out the numerical analysis of the SR one usually performs the Borcl transformation on both sides 
of the sum rule 



bu(q 2 )^ lim ^—f-{Q 2 ) n+1 [ -77x2) n(<3 2 ) 



Q2/n=M 2 =fixed 



where M is called the Bore l ma ss which represents the scale r = 1/M in Euclidean time between the two currents 
in the correlator, Eq. (|4"3"1) .[145] The Borel transformation of the sum rules has two major advantages. The first, 
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it produces an enhancement of the glueball pole contribution in the right hand side of the SR and, the second, it 
suppresses, in the theoretical left hand side of the SR, the contribution from the higher power corrections, which arise 
from high dimension condensates. 

The main point of the QCD SR philosophy is the implementation of the interaction of the high virtual valence 
quark-gluon system with the soft vacuum quark and gluon fields, whose strength is determined by the values of the 
vacuum condensates. The condensates carry very important information about the long-range quark and gluon field 
fluctuations in the QCD vacuum. The usual assumption is that one can calculate the hard part within perturbative 
QCD and the soft part can be parameterized in terms of condensates. It has been demonstrated that, within the 
OPE approach, in general, the contribution of only a few low dimension condensates is sufficient to describe rather 
well the properties of the ground hadronic states. [72J However, it was found, many years ago, that in some specific 
channels, which include the quark-gluon subsystem in a spin zero state, the standard OPE expansion does not work 
and one needs to incorporate more precise information about the structure of the short-range gluonic fluctuations 
in the QCD vacuum. [63, [zH [zl| The most promising candidate for those short-range fluctuations are instantons 
(see ref. (z3|). The instantons describe tunneling processes which rearrange the QCD vacuum topology in localized 
space-time regions small enough to affect the x-dependence of the correlators over distances |x| <C Aqq D . The 
average size of instantons in the QCD vacuum is small p c « 0.3 fm. Therefore, the energy scale related to instanton 
effects is rather large, raj — l/p c ~ 600 MeVw 3Aqcd- Due to such large scale, instantons influence strongly not 
only the dynamics of the low mass pseudoscalar and scalar mesons, but also the properties of the baryon octet and 
multiquark states, ffdi, l77l. |78| It has been shown that specific instanton induced quark-quark, quark-gluon and gluon- 
gluon interactions arising from intermediate distances between hadron constituents can be responsible for various 
observed features in hadron spectroscopy and hadron reactions. [7^, [7^, HH The new ingredient within QCD SRs 
related to instantons, the so-called direct instanton contribution, provides exponential terms, ~ exp(— 2Qp c ) in the 
expansion of the correlator, in addition to the power terms, ~ \/Q n , arising from the standard OPE. 




(a) (b) (c) (d) (e) 



FIG. 6: The diagrams (a) and (b) represent pQCD contributions, diagram (c) represents the contribution arising from the gluon 
condensate, diagram (d) is the contribution from three-gluon condensate and diagram (e) is the direct instanton contribution. 

Finally, the theoretical part of the SR includes the following terms (see Fig. [6]) 

n thoor (Q 2 ) = n port (Q 2 ) + n cond (Q 2 ) + n inst (Q 2 ). (50) 

The first two terms in Eq. (|50l) are calculated using Feynman rules giving for the perturbative part for the correlator 
of the gluonic field strength [7l| 

(TG« v (x)G» af} (0)) peit = -i5 ab J J^r^CpJe-**, (51) 

where 

IV Q /3(p) = PuPuQaP + PvPpg^a ~ P^Vjidva ~ PuP a g^f3- 

The term associated to the vacuum condensates is evaluated in the so-called fixed point gauge. In this gauge the soft 
vacuum gluonic field can be represented by the field strength 

K( x ) a ~ %7 G lu M- ( 52 ) 

M v 'cond 2 ^cond v ' v ' 

Using that the QCD vacuum is colorless and Lorenz invariant, and Eq. (|52p . one obtains the contribution of the 
vacuum condensates to the OPE for the gluonic channels. It turns out proportional to the value of gluon condensates, 

(a s G 2 ) = (a s G^G« u ), (g s G 3 ) = { 9s f abc G%Gl p G c ^) . (53) 

Unfortunately, the values of the gluonic condensates are not well known. Even the lowest dimension gluon condensate 
(a s G 2 ) w 0.035-0.075 GeV 4 is fixed by considering the SRs for quark systems and has large uncertainties (see discussion 
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in ref. [82J ). Other gluonic condensates have been estimated using models for the QCD vacuum. For example, the 
tree gluon condensate within single-instanton approximation is quite small (g s G 3 ) « 0.27 GeV 2 (a s G 2 ). Lattice 
calculations are also uncertain. For the lowest dimension condensate the quenched value, 0.14 ± .02 GeV 4 , is about 
one order of magnitude bigger than the unquenched one, 0.022 ± .005 GeV 4 . [HI The uncertainties in the condensate 
values lead to effects of the order of at least 20% in the extracted masses for glueballs. 

The direct instanton contribution is calculated by going to the Euclidian space-time and by introducing in the 
correlator (|4"5)) the instanton field strength, which in the regular gauge has the following form 

2 

u^yx, x a — — -75— — ~y>, v 04 / 

g s ((x - xq) z + p z y 

where xq is the position of center of instanton and rja^ are the numerical t'Hooft symbols. To get the final result one 
should also integrate over instanton position, size and orientation. 

Let us study in some detail the scalar glueba ll g roundstate, J = ++ , which would be an ideal state to find 
since it has many fundamental connotations. [661 184. l85l I86L l87j In gluodynamics, as we have seen, the situation that 
arises from lattice calculations is clear and the masses of the scalar glueballs are large m > 1 GeV. However, when sea 
quarks are considered no firm conclusion about the scalar spectrum can be drawn. The theoretical calculations based 
on QCD SRs and/or low energy theorems lead to contradictory results. Its properties, i.e., mass, decay channels 
and widths still differ among the various approaches. Dominguez and Paver [88| . Bordes, Peharrocha and Gimenez 
[89| . and Kisslinger and Johnson [9(|, obtain, using low energy theorems and/or SR calculations with (or without) 
instanton contributions, a low lying and narrow scalar glueball (mass < 700 MeV and < 100 MeV). Narison 
and collaborators [U using two (substracted and unsubstracted) SRs, get a broader (200-800 MeV), heavier 
(700-1000 MeV) scalar glueball in this range, whose properties imply a strong violation of the Okubo-Zweig-Iizuka's 
rule. In a recent SR calculation, Forkel [93 |. obtains the scalar glueball at 1250 ± 200 MeV. However, he has some 
strength at lower masses which he is not able to ascribe to a resonance in the fits. 

Having made these comments let us now see how a calculation proceeds and we let the details of the various 
approaches just described for the reader to make his own opinion. The calculation we describe is in line with the work 
of Forkel. 

In this case the perturbative contribution arising from Fig. [SJa) and (b) is proportional to a high power of the Q 

n pcrt (Q2) = ( ^ )2Q 4 log( Q2 /M 2 )(Ai + (55) 

where Aj are some numbers. Therefore, one needs to apply three substractions to have convergence for the dispersion 
relation 94 



n*w (Q2) = n(0) + Q2n , (0) + V n „ (0) _ QP_ r ds Im "^) ( 56) 

The SRs for the following Borel transforms 

B k {t) = -B[(-l) K Q 2K n(Q 2 )] (57) 



arc 



-1 />00 7 

B-x = -n(0) + - / -e- ST Imn(s) 

T Jo S 



B K = — I dss k e^ ST lmU(s), k > -1. (58) 



After moving the contribution from continuum to the left hand side (theoretical part) of the SR we have 

1 f 3 ° ds 

sr_i(r,s ) = -n(0) + - / - e - ST imn( s ) 

x Jo s 



Finally the SRs become 



5 K (r,s ) = -/ dss k e~ ST lmn(s), k>-1. (59) 



S-i(t,s ) + H(0) = ^/S^G,e mG . T 

i 

5 K (r, So ) = ^2f^m 2 G T 4 e m ^ T , k > -1. (60) 
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Usually the lowest n value SRs are used to extract the value of glueball mass because they are more stable with 
respect to the continuum threshold 



m 2 



Si(t, s ) 
Sq(t, s ) 



2 Sq(t,S ) 

= s,(r,„) + n(o)' (61) 

In principle, one should obtain the same mass for the glueball in the above two SRs. It was found that the extracted 
masses for the scalar glueball are quite different in the two. However, as a recent analysis shows, if one includes the 
direct instanton contribution the mass arising from the two SRs is the same. [lH [9j| [96j 

These works also suggest a prominent role of the instantons in the binding of the ++ glueball and show relations 
between the main properties of the + + glueball (mass and size), and the bulk features of the instanton distribution. 

In the meantime, the study of instanton contributions to SRs has been extended to other glueball channels. [93|, [§6| 
While direct instanton contributions are expected to be small in the tensor channel (mainly due to the absence of the 
leading correction with respect to the small packing fraction of instantons in the QCD vacuum), they turn out to be also 
important in the pscudoscalar channel. However, we should stress that the instanton effects are counterbalanced in the 
pseudoscalar channel by a screening of the topological charge. [931] This screening occurs when the quark contribution 
is important and leads to a prediction consistent with a low energy theorem. The most accurate result , along these 
lines, for the masses of the scalar (S) and pseudoscalar (P) glueballs including direct instanton effects in traditional 
SR calculations is 

m s ~ 1-25 GeV, m P » 2.2 GeV. 

Within the uncertainties of the QCD sum rule approach these values are in agreement with lattice calculation, which 
are quenched calculations. Note that the SR predictions include sea quark effect and thus one should be very careful 
when comparing with the quenched lattice results. In the quenched approximation, the topological charge screening 
disappears and the instanton effect in the SR for the pseudoscalar glueballs is bigger. 

We should emphasize that one central problem of glueball spectroscopy, namely the mixing between quarkonium 
and the spin zero glueball states still has not been solved so far. Recently, some steps have been done in t his 
direction, [U [9§] and it has been shown that the instantons induce a strong mixing between the two states. [l46j ] In 
this respect we would like to point out that from our point of view, the tensor 2 ++ glueball channel is the most clean 
gluonic channel. Indeed, for this state the leading instanton contribution is zero, as follows from the structure of the 
interpolating current for this channel, which is proportional to the energy- momentum tensor, Eq. (|46[) . and from the 
fact that the instanton is a vacuum gluon field with zero energy. Therefore, one can expect very small tensor glueball 
mixing with quarkonium and, as a consequence, a small width for this state and thus the possibility t o separat e it 
from the tensor quarkonium states. The SR prediction for the mass and decay constant of 2 ++ state is [lOOL Il0l| 

m T w 2.00 GeV, T t ^ + kk+ wi < (119 ± 36) MeV. 

Due to absence of big uncertainties in the direct instanton contribution and the small expected mixing with quarko- 
nium states in this channel, this SR prediction is on more solid ground than the those for the zero spin glueball 
states. 



V. GLUEBALL PRODUCTION AND DECAYS 



One of the most striking features of QCD is the prediction that glueballs might exist, a prediction which has proven 
difficult to verify. We have seen how different theoretical approaches describe their properties, and now we are going 
to discuss possible scenarios for their experimental production and detection. 

The strategy relies on a very few assumptions. Glueballs are extra states, beyond the qq spectrum. To exploit 
this we must understand the "ordinary" qq spectrum very well, using data from J/ip, B, and Z decays, and from 
pp, irp, 77, and 7^ scattering. Glueballs are flavor singlets so their decays should be SU(3)f symmetric, if the 
effects of violation of SU(3)f symmetry are small. Since hybrid and qqqq states are also "extra" states and some are 
flavor singlets, the only distinguishing property unique to glueballs is their strong coupling to the color singlet digluon 
channel. 

A convincing way to see a glueball would be to detect a resonance whose quantum numbers are not possible for 
mesons composed of quark and antiquark, e.g. J PC = 0~~,0 H ,1 h ,2 H . These glu eball s are sometimes called 
oddballs. They are predicted to be rather narrow and easy to identify experimentally. [102J These states seem to 
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appear for very large masses. The lightest oddball wit h J P C — 2 H and a predicted mass of 4.3 GeV will be within 
range of the future experimental program at PANDA. [l03j ] It is conceivable that comparing oddball properties with 
those of non-exotic glueballs will reveal deep insight into the glueball structure since their spin structures are different. 
However, if oddballs happen to appear at lower masses we cannot rule out a glueball interpretation but they should 
be strongly mixed with hybrid mesons. 

However, even non exotic glueballs can be identified by measuring an overpopulation of the experimental meson 
spectrum and by comparing masses, quantum numbers and decay channels, with predictions from models or lattice 
QCD. The best scenario is to look for Zweig forbidden processes, since the decay into glueballs is dominant in these 
cases. |l04j These processes are known under the name of gluon-rich processes, some are depicted in Fig. [7J Let 
us specify some of the processes that fulfil these requirements and are optimal experimental scenarios to search for 
glueballs: 

(i) J/ip decays: the most suggestive process is the radiative J/ip decay. The J/ip is narrow; the DD threshold is 
above the mass of the J/ip and the OZI rule suppresses decays of the cc system into light quarks. In most decays, 
the J/ip undergoes a transition into 3 gluons which then convert into hadrons. This is a high multiplicity process 
difficult to detect. But J/ip can also decay into 2 gluons and a photon. The photon can be detected, the two 
gluons interact and must form glueballs. BES III will provide huge J/ip data samples allowing definitive s tudi es 
of 3/tp decay and, especially, partial wave analysis of the glueball-preferred radiative J/ip decay channel. jl05l | 

(ii) Central production: glueballs decay into hadrons and hence hadro-production of glueb alls is always possible. 
Central production is a process in which glueballs should be abundantly produced. [106j In central production 
two hadrons scatter diffractively, gluon-gluon fusion processes should be abundant and no valence quarks are 
exchanged. The process is often called double pomeron exchange, with the pomeron presumed to be a multi- 
gluon color singlet. The absence of valence quarks in the production process makes central production a good 
place to search for glueballs. 

(iii) pp annihilation: quark-antiquark pairs annihilate into gluons, they interact and may f orm glueballs, e.g. 
pp — > 7r°G. A new era of precise experiments will start when PANDA starts to operate. 103]. The energy range 
of this experimental program will allow the study of oddballs. Glueballs can either be formed directly in the 
pp-annihilation process, or produced together with another particle. In both cases the glueball decay into final 
states like (pep or cprj would be a favorable reaction below 3.6 GeV, while J/iprj and J/ip(p are the first choice for 
more massive states. 

(iv) 77 fusion: production of glueballs should be suppressed in this case since photons couple to the intrinsic 
charges. So we should expect a glueball to be strongly produced in radiative J/ip decays but not in 77 fusion. 
Radial excitations might be visible only weakly in J/ip decays but they should couple to 77. However, if a very 
low mass glueball would exist, this would be one of the cleanest ways to find it in a difficult experiment due to 
the low counting rate. 

(v) photoproduction: it has been a continuous source of meson resonances. Photoproduction is expected to 
be particularly effective in producing exotic hybrids and glueballs. The photon, via Pomeron exchange from 
the vacuum, or via its hadronic component, can create states with exotic J PC . These mesons can be hybrids, 
glueballs or mixed states. Moreover, also non exotic hybrid (glueball) states can be produced which manifest 
themselves as extraneous states that cannot be accommodated within qq nonets. However, there is little data on 
photoproduction of light mesons. This will drastically change when the Glu eX experiment goes into operation, 
since it is designed to collect data of unprecedent statistics and quality. jl07l | Moreover, the (linear) polarization 
of the beam will allow the identification of (exotic) quantum numbers and the determination of the details of 
the production mechanism of (glue) mesons. 



The status of glueball observation would become clearer with a combination of data of pp and e + e~ machines with 
large statistics, however since the glueball width into two-photons is small, this would require very high luminosity. 
It has been pointed out that heavy-ion colliders, due to their large center of mass energy, allow access to the Regge 
region s 3> \t\ and hence to the p roduction of g lueballs in peripheral collisions through photon-photon collisions as 



Distinctive features can be derived from the decays of the glueballs since they are flavour singlets: decays to r\Tj 
identify a flavour octet and radiative decays of glueballs are forbidden. However, these arguments have to be taken 
with care since mixing of a glueball with mesons, having the same quantum numbers, can occur and would dilute any 
selection rule. 
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FIG. 7: Diagrams possibly leading to the formation of glueballs: radiative J/ip decays, pomeron-pomeron collisions in hadron- 
hadron central production, and in pp annihilation. 



An interesting idea recently proposed is that of chiral suppression. [1 1 1| If chiral symmetry breaking in glueball 
decay is dominated by quark masses, then the coupling of a spin zero glueball to light qq pairs is chirally suppressed. 
An interesting consequence is that radiative J ftp decay becomes a filter for new physics in the J — channel, since 
at leading order radiative decays to spin zero light quark mesons are suppressed while radiativ e decays to J = 
glueballs, hybrid, and four quark states are not suppressed. This idea is not without controversy. [ll2l Ill3j | 

The glueball groundstate J PC = ++ would be an ideal state to find since it has many fundamental connotations. [6tj 
l&H . [85l. l86l. H3] In gluodynamics, as we have seen, the situation that arises from lattice calculations is clear and the 
masses of the scalar glueballs are large m > 1 GeV. However, when sea quarks are considered no firm conclusion 
about the scalar spectrum can be drawn. The theoretical calculations based on QCD sum rules and/or low energy 
theorems lead to contradictory results. Its properties, i.e., mass, decay channels and widths still differ among the 
various approaches. 

Experimentally there are too many isoscalar, scalar mesons between 0.6 and 1.75 GeV to be explained by the naive 
quark model alone. From the theoretical point of view there has been a lot of debate about the structure of these 
states. Let us summarize here the different theoretical e xplan ations for the spectrum. Jaffe obtains the /o(600) and 
the /o(980) in the bag model as cryptoexotic qqqq states. [114| The p-wave qq scalar nonet is likely to lie in the region 
of the other spin-triplet p-wave nonets, with isoscalars roughly between ~ 1250 and ~ 1600 MeV. [l3[ Dominguez and 
Paver. [88) Bordes, Peharrocha and Gimenez.[89j and Kisslinger and Johnson. [90| obtain, using low energy theorems 
and/or sum rule calculations, a low lying and narrow glueball (mass < 700 MeV and r w7r < 100 MeV). Vento using low 
energy theorems has also proposed the existence of a low mass (m < 700 MeV) and almost sterile scalar glueball. [87| 
Narison and collaborators. [9ll 92 1 using sum rules, get a broader (200-800 MeV), heavier (700-1000 MeV) scalar 
glueball in this range. Forkel, 93[ obtains a scalar glueball at 1250 ± 200 MeV with a large width (^300 MeV). 
Between ~ 1400 and 1750 MeV there are three I,J PC = 0,0++ states: / (137Q), / (1500), and / (1710). In the 
analysis of Vento [§7| two scalar glueballs appear in this range, an intermediate mass glueball (~ 1300 MeV) and 
another one in the upper range (~ 1500 — 1700 MeV). The narrow state at 1500 MeV, discovered in pp an nihilation 
by Crystal Barrel (CB) is considered by Amsler and Close a good can dida te for the glueball groundstate [l3l ITTsl 1 1 1 6j 
while Chanowitz considers the /o(1710) as the scalar groundstate (lil fTTl| . Finally, that the scalar glueball is shared 
by both resonances is also contemplated. jl4] 



VI. GLUEBALLS AND THE QUARK-GLUON PLASMA 

Relativistic heavy-ion collisions might be a tool to produce glueballs. It is conceivable that glueball production 
becomes a dominant part in central nucleon collisions. 
Two scenarios will be analyzed: 

(i) Quark Gluon Plasma phase: one expects that at some point ab ove a certain critical temperature a plasma 
of quarks and gluons, named quark-gluon plasma (QGP) is formed. [ll7| | This phase is characterized by a large 
amount of thermal gluons. As the QGP cools down the gluons can form singlet configurations via the color 
interaction. 

(ii) Strong Coulomb phase: a recent formulation of the dynamics in the region above the transition temperature 
T<7, based on a description of recent experiments in ultra relativistic heavy ion collisions. [118| | states that, despite 
de-confinement, the color Coulomb interaction between the constitue nts is strong and a large number of binary 
(even color) bound states, with a specific mass pattern, are formed. ( 1 1 Oj ] This phase we call Strong Coulomb 
Phase (SCP). The QGP phase occurs at a much higher temperature Tqqp > (2 — 3)Tc when the bound states 
dissolve. 
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Let us describe some attempts to find signatures arising from QGP. The basic idea is that during the hadronization 
process in this gluon rich environment the gluons combine due to the strong force into glueballs. As the QGP cools 
further and transforms into hadronic matter, these glueballs decays into conventional hadrons giving rise to signatures 
of their existence. It has been claimed that the existence of glueballs alters the K/n ratio in the final state. [120J 

It is conceivable that glueball production becomes dominant in central nuclear collisions since the existence of QGP 
provides a gluon rich environment especially at high energy density. In this scenario the lowest mass glu eball s are 
copiously produced. Particular decay modes ++ , 2++ — > KK and ++ — > ir + n~l + l~ have been investigated jl2l| and 
enhancements associated with possible glueball production observed. Search strategies, including dipion production 
have been also proposed. [l22j | 

Estimates within thermal models for the multiplicities of scalar glueballs in central A u-Au collisions at present 
and future experimentally available energies, i.e. from AGS to LHC, have been presented. |l23l ] For the experimental 
identification of glueballs one can use the decay modes G — ► KK, G — ► 77 and G — ► 2irl + l~ . Despite of small 
branching ratio, the 27 channel has the important advantage that photons have practically no rescattering in the 
hadronic medium. This analysis leads to maximal multiplicities for the glueball of 1.5 — 5% of the <t> mes on multiplicity. 
Even larger yields are expected in the case of explosive hadronization of the quark gluon plasma. |l23l | 

Let us now turn to the second sce nario . From its inception two decades ago, the high-T phase of QCD commonly 
known as the Quark-Gluon Plasma, 117j | was described as a weakly interacting gas of "quasiparticles" (quarks and 
gluons). Indeed, at very high temperatures asymptotic freedom causes the electric coupling to be small and the 
QGP to be weakly interacting. At intermediate temperatures of few times the critical temperature T c of immediate 
relevance to current experiments, there is new and growing evidence that the QGP is not weakly coupled. In this 
region QCD seems to be close to a strongly coupled Coulomb regime, with an effective coupling constant a f=a 0.5-1 
and multiple bound states of quasiparticles. [TT^ | This phase we call Strong Coulomb Phase. This de scription is not 
universally accepted since some lattice calculations do not to find bound st ates a bove the transition. jl24 H25| Also, 
one must be aware that other explanations of the data have been presented. [l26| In order to clarify these issues it is 
interesting to find clear experimental observables that would shed some light into the discussion. 
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FIG. 8: The solid (dotted) line is the temperature dependence of scalar (pseudoscalar) glueball mass (left). The temperature 
dependence of the gluon condensate at T > T c . The solid line is total condensate, the dashed (dotted) line is the scalar 
(pseudoscalar) glueball contributions (right). 



Kochelev and Min consider the properties of scalar and pseudoscalar glueballs in SCP. [l27j | In an effective La- 
grangian approach, based on the low-energy QCD theorems, they find out that scalar glueball remains massive above 
deconfinement temperature (see Fig. [5]). At the same time, pseudoscalar glueball changes its properties in QGP 
in a drastic way. Indeed, this glueball becomes massless at T > T c and therefore it can contribute strongly to,the 
bulk properties of the plasma (see Fig. [8|). They demonstrate that the disappearance of pseudoscalar glueball mass 
above the deconfinement temperature and its strong coupling to gluons gives the rise to the sign change of the gluon 
condensate in the pure SU(3) C gauge theory as observed in the lattice calculations at T ps T c (see Fig. I5]). |l28l | The 
strong nonperturbative coupling of the glueball to the gluons leads to the conjecture that one might expect that the 
role of very light pseudoscalar glueball in QGP must be quite similar to the role played by the massless pion in nuclear 
matter below deconfinement temperature. 

In SCP despite de-confinement the color Coulomb interaction between the constituents is strong and a large number 
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of binary (even color) bound states, with a specific mass pattern, are formed. [11 Of ] Wi th th is input, the scenario 
envisaged by Vento for gluodynamics, i.e. the theory with only gluons goes as follows. [l 2 9j | The strong Coulomb 
phase is crowded with gluon bound states and the lightest is the scalar glueball, labelled g. As one moves towards 
the dilution limit, the binding energy of these states decreases, the gluon mass increases, and t heref ore the color and 
singlet bound states increase their mass softly until the gluons are liberated forming a liquid. [119( However, as the 
system cools towards the confining phase, color and singlet states decay into the conventional low lying glueballs, in 
particular g. Thus the number of g's becomes large. 

This reasoning generalizes to QCD since in the SCP the multiplicity of glueball channels is larger than in the 
confined phase. The ratio of glueball to meson channels goes from 1 to 8 below the phase transition to 1 to 2 
above. [TT!| Thus the number of scalar glueballs is much larger in SCP than in the cold world. 

As the fireball cools a "large number" of gluonic bound states decay by gluon emission into <?'s. The emitted gluons 
form new bound states of lower mass due to the strong color Coulomb interaction. As we approach the confinement 
region the mass of the color bound states increases and it pays off to make multiparticle color singlet states, which 
decay by rearrangement into ordinary color singlet states. Since the coupling is strong and the phase space is large, 
these processes take place rapidly. Thus in no time, close to the phase transition temperature Tc, a large number of 
scalar glueballs populate the hadronic liquid. 
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FIG. 9: Behavior of the masses of a,n and g across the QGP phase transition according to model calculations. g± label the 
upper(lower) limits of the g in mass model calculations. 130 



Let us now analyze the experimental signal of these phenomenon. If we assume that the a is the 0(4) partner of 
the 7T in the chiral symmetry realization of QCD, its mass decreases when approaching the phase transition, becoming 
degenerate with the pion at Tc (see Fig- EJ) - Beyond Tc, in the SCP, chiral symmetry is restored, and tt and a remain 
degenerate for T > Tc- Thus in the SCP the a can only decay in 27 for obvious kinematical reasons. The glueball 
g does not vary its mass in this region appreciably. Thus even before we reach Tc, the mixing between g and a 
disappears (see Fig. [9]) and g becomes stable around Tc- However, in the SCP the mass of the a increases and in a 
certain region of T it again becomes degenerate with g and mixing is restored. Thus the physical g is able to decay, 
once the a component is attained, only to 27. 

The enhancement in the number of gs with respect to the hadronic phase arises because of the larger population of 
glueballs in the SCP, as described above, and because these particles are stable in the medium against the dominating 
hadronic decays (see Fig. [9|) . Thus a clear signal for the existence of a low mass scalar glueball and a confirmation of 
the SCP scenario would be two 27 peaks corresponding to g and the er-meson as shown in Fig. [TO] 



The investigation on jets in th e relativistic heavy ion collisions at RHIC provides a deep insight into the properties 
of the quark-gluon plasma. jll8| One of the important RHIC discoveries is the jet quenching phenomenon coming 
from the partonic energy loss in QGP. In the conventional approach to the jet quenching the perturbative (pQCD) 
type of energy loss is taken into account by the channels, elastic and radiative, of one-gluon exchange between the jet 
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FIG. 10: Expected fit to the two-photon invariant mass spectrum in central Pb-Pb collisions after substraction from the 
background. The 27 decays should allow for a clear separation of g and a. The figure includes an estimate of the effect of the 
hadronic widths on the fits (dot-dashed line). 



and the massless gluons and quarks, fl 3 1| However, the large quark-gluon rapidity density dN qg /dy s» 2000 which is 
needed to describe the RHIC jet quenching data within this approach, seems to be in contradiction with the restriction 
dN q g/dy < 1/AdS/dy « 1300 coming from the measured final entropy density dS/dy w 5000. 132] Furthermor e, th e 
lattice calculations show that even at very high temperature gluons and quarks still interact strongly in QGP. [1281 ] 
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FIG. 11: The temperature dependence of gluon energy loss (left). The temperature dependence of quark energy loss (right). The 
solid (dashed) line is glueball (pQCD) contribution . 



As seen above, it has been suggested that the glueballs, the bound states of gluons, can exist abov e deconfinement 
temperature and may play an important role in the dynamic of strongly interacting QGP. [1271 Il29j ] In particular, it 
has been suggested that a very lig ht ps eudoscalar glueball can exist in QGP and might be responsible for th e res idual 
strong interaction between gluons. [l27j ] The lattice results showing a change of sign of the gluon condensate [l28| and 
a small value of the topological susceptibility above T c can be explained in the glueball picture as well. Furthermore, 
one expects that the suppression of the mixing between glueballs and quarkonium states in the QGP leads to a smaller 
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width for the former as compared to the vacuum. [l29| | This property opens the possibility for a clear separation of 
the glueball and the quark states in heavy ion collisions. Such separation is rather difficult in other hadron reactions 
due to existence of strong glueball-quarkonium mixing in the vacuum. 

Min and Kochelev made an estimate of the energy loss induced by interaction of an energetic parton, whic h was 
produced in the hard scattering of two heavy ion's partons, with glueballs in the hot quark-gluon plasma. [1331 ] They 
showed that such contribution leads to a significant energy loss as can be seen from Fig. [11] In particular, for th e glu on 
jet such contribution is about a few GeV/fm and approximately twice larger than the perturbative elastic loss. [134j It 
should be pointed out that more than one half of contribution to the gluon energy loss comes from interaction of gluon 
with the light pseudoscalar glueball in QGP. Thus in spite of the fact that for the quark jet the glueball contribution 
is smaller than perturbative elastic loss, it can not be neglected in comparison with latter one. The existence of such 
light bound state of gluons above T c is crucial for the understanding of the large observed partonic energy loss in 
QGP. Therefore, not only pQCD type of energy loss but also glueball-induced loss, arising from existence of scalar and 
pseudoscalar glueballs in QGP, are important for the understanding of the RHIC results such as the jet quenching. 

VII. OTHER DEVELOPMENTS 

A. Glueballs and the Pomeron 

It is well known that in the many high energy reactions with small momentum transfer the exchange by the 
highest-lying Regge trajectory, called as a soft Pomeron, gives the dominant contribution. This exchange carries 
vacuum quantum numbers and has very peculiar properties in comparison with the usual Regge pole trajectories as 
the p, 7r, etc. From the analysis of the varions cross sections it follows that the Pomeron trajectory has a linear 
behavior 

J(t = to 2 ) = 1.08 + 0.25m 2 . (62) 

The Pomeron does not seem to be related to the usual mesons since the latter have usually lower intercepts and very 
different slopes. There has been a long-standing speculation that the physical particles on the trajectory might be 
glueballs. 

Meyer and Teper investigated the pure gauge spectrum on a lattice to check its compatibility with a linear 
trajectory. [l8j They calculated the lightest states J PC = 2 ++ ,4 ++ ,6 ++ since the trajectory has an even signa- 
ture. Their masses (within errors bars) are consistent with the Pomeron trajectory (|62p. These studies of quarkless 
QCD helps to understand the close relation between the Pomeron and the pure glue states, but it is worth mentioning 
that the physical particles should probably mix with quarkonia. The linearity of the trajectory, supported by these 
arguments, favors flux tube and stringy pictures as discussed in the constituent model section, at least for (even) ++ 
glueballs. 

B. Glueball - qq mixing 

Most of our review has been dedicated to the study of pure glueball states. We have assumed in our discussion that 
glueballs arc physical objects, which is fine for the purpose of the developments thus far. But in order to establish 
contact with reality we have to discuss in some detail mixing, an idea which has been mentioned in passing on several 
occasions. We first discuss the mathematics of mixing in a two state model and then apply the reasoning to the 
physical reality in the isosinglet scalar sector. 

Mixing arises when the total Hamiltonian is not diagonal in the glueball and qq Fock space. In this case the physical 
states come out as superpositions of glueballs and qq states. Let us assume that we have a pure glueball, we label 
g, and a qq state, we label a, which have the same quantum numbers and, for simplicity, also the same mass to, in 
a certain approximation. We next relax this approximation and the Hamiltonian in the reduced g and a Fock space 
becomes 

( 5 to + Ato ) ( 63 ) 
The diagonal basis of this Hamiltonian can be presented as, 



\g) = cos (0/2) |<?> -sin (0/2) |<t), 
[or) = sin(0/2)|.g) + cos (0/2)|cr), 



(64) 
(65) 
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FIG. 12: The limiting values for the masses of the physical g (solid-dashed lines) and a (solid-short-dashed lines) are shown 
as a function of the mixing angle. 



where the tilde labels the physical particles and 9 is the mixing angle. 
The masses of the physical particles become 



where 



and 



Am Am 
rrig = m H r, ms = m H — h r, (66) 



9 ^ 

tan6>=— - (67) 



Am 




Am .„„. 
= 2^0- (68) 

In Fig. [T^Jwe represent the masses of the physical states as a function of the mixing angle 9. The curves separate 
possible mass regions. The curves show that the two state mixing scenario for positive Am leads to a "light" glueball 
and heavier meson. The opposite result can be obtained with a negative Am. The procedure can be generalized to 
more states and to different initial masses as we next see. This mechanism has been used to try to disentangle the 
complex isoscalar sector. 

Vento has shown that a low mass, almost invisible isosinglet scalar glueball, could hide in the low energy spectrum 
due to small mixing with the /o(600), and pointed out the possible existence of other glueballs components among 
the remaining isosinglet scalars /o(1370), /o(1500) and fn(1710).[87j 

Close et al. have discussed the content of the the isosinglet scala r mesons ,fo(1370), /q(1500) and /o(1710) in 
an attempt to discover which of them is dominantly a scalar glueball. [135l Il36| These authors have suggested that 
/o(1500) is primarily a scalar glueball, due partly to the fact that /o(1500), dis cove red in pp annihilation at LEAR, 
has decays to 7777 an d nn' which are relatively large compared to that of 7T7r [l37j | and that the earlier quenched 
lattice calculations d, Il38j | predict the scalar glueball mass to be around 1550 McV. Furthermore, because of the 
small production of itit in /q(1710) decay compared to that of KK, they claim that the /o(1710) is primarily ss 
dominated. In contrast, the smaller production rate of KK relative to irir in /o(1370) decay leads to the conjecture 
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that /o(1370) is governed by the non-strange light quark content. Based on these observations, they have proposed 
a flavor- mixin g scheme to consider the glueball and qq mixing in the neutral scalar mesons /o(1710), /o(1500) and 
fn(1370). |l35f Fits to the measur ed scalar meson masses and their branching ratios of strong decays have been 
performed in several jl35l 1 136L Il39j | leading to a mixing matrix of the form 

/o(1370) \ / -0.91 -0.07 +0.40 
/o(1500) = -0.41 +0.35 -0.84 
/o(1710) J \ +0.09 +0.93 +0.36 

where \N) and 15) denote the quarkonium states (\uu) + \dd))y/2 and \ss), and \G) denotes the pure scalar glueball 
state. Thus for Close et al. the /o(1500) is composed primarily of a scalar glueball. 

Another analysis has been carried out by Cheng et al.. 140] Two lattice results are employed as the starting point; 
one is the ap prox imate SU(3) symmetry in the scalar sector above 1 GeV for the connected insertion part without qq 
annihilation, [14l| and the other is the scalar glueball mass at 1710 MeV in the quenched approximation. [I?], [H[ In the 
SU(3) symmetry limit, /o(1500) becomes a pure SU(3) octet and is degenerate with ao(1450), while /o(1370) is mainly 
an SU(3) singlet with a slight mixing with the scalar glueball which is the primary component of /o(1710). These 
features remain essentially unchanged even when SU(3) breaking is taken into account. The observed enhancement 
of w/o(1710) production over <^>/o(1710) in hadronic J/ip decays and the copious /o(1710) production in radi ative JJjb 
decays lend further support to the prominent glueball nature of /o(1710). Furthermore, chiral suppression [H [m, E3 
is advocated to obtain the following mixing matrix, 

' / (1370) \ / +0.78 +0.51 -0.36 
/ (1500) = -0.54 +0.84 +0.03 
/o(1710) J \ +0.32 +0.18 +0.93 

Therefore for Cheng et al. it is the /o(1710) is the particle composed mostly of a glueball state. Thus a definitive 
conclusion on this problem is still lacking. 

Finally we would like to mention that the quarkonium-glueball mixing gives also a strong influence to the properties 
of the lowest mass pseudoscalar glueball. [97| Il43j 




VIII. CONCLUSIONS 



The leitmotif of this review has been that the study of glueballs, states where the gauge field plays an important 
dynamical role, is important to understand the nonperturbative behavior of QCD. This study requires precise and 
abundant experimental input, which despite a lot of effort has not achieved a level of understanding that allows a 
unique theoretical interpretation. The lack of clarity arises from the fact that the theoretical developments are not able 
to determine in a well defined manner how the production and decay properties of glueballs are distinct from those of 
conventional mesons. But the future is bright for several reasons. New energy domains will be studied by BESIII and 
in the future by PANDA, and therefore the possibility of producing oddballs, glueballs with exotic quantum numbers 
will open. Moreover, at low energy, both Crystal Barrel and in the future GlueX will produce light mesons with a 
level of statistics, and precision for analysis, that will be able to separate exotic from non exotic behaviors and quark 
model nonets from particles extraneous to them. In this process the new theoretical developments will be crucial to 
guide and/or interpret the results. 

In this review we have discussed three major approaches, lattice QCD, QCD based constituent models and QCD 
sum rules, with a short excursion into the AdS/QCD formalism, which at present leads ultimately to a constituent 
model type of description. 

Due to the lack of observable states, the lattice QCD results have played the role of experimental data. Lattice 
QCD is a powerful technique, especially for the determination of masses. The spectrum of the pure gluodynamics (the 
pure gauge theory), equivalent to the so called quenched approximation of lattice QCD, is well established and the 
calculational errors are under control. However the various treatments of quark loops (unquenched lattice QCD) and 
the effects of the mixing of glueballs with mesons are still a matter of debate and no firm conclusions can be drawn 
for the real world. Moreover, important features of the study would be the determination of the decay properties of 
the various glueball candidates, and lattice QCD being a theory described in Euclidean space-time does not possess 
the capability of describing asymptotic states and therefore of determining decay properties. However, a full lattice 
QCD calculation with unlimited and precise calculational power should produce the experimental spectrum, if as we 
strongly believe QCD is the theory of the hadronic interactions, and moreover, if one can extend lattice QCD into 
Minkowski space-time one would be able to determine the decay properties. 
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Given the lack of experimental and detailed theoretical knowledge, models turn out to be an interesting laboratory 
to test ideas of the various perturbative and nonperturbative mechanisms in the theory. We have discussed in some 
detail constituent gluon models for glueballs from different perspectives. Most of the pioneer work considered glucballs 
as strongly bound states of two (or three) heavy spin-1 gluons. The comparison with lattice QCD shows that the 
spectrum obtained by these approaches does not correspond to lattice QCD. However, if one implements a formalism 
dealing with only transverse gluons, the spurious states (induced by the gluonic longitudinal components) disappear 
and the hierarchy in the spectrum of lattice QCD is recovered. These improvements have been developed up to now 
only in the two-gluon sectors and the generalization to negative C-parity is still lacking. It remains for the future to 
find out which constituent glucball model confirms the lattice QCD spectrum for the higher states. 

The dynamics in these models is described, in general, by a linear plus Coulomb potentials. This dynamics, however, 
does not remove the degeneracy of the scalar and pseudoscalar glueballs. The candidate instanton induced interaction, 
which arises naturally in QCD and which has very particular properties under spin-parity, is the natural candidate to 
do the job. In the quenched approximation, the instanton induced interaction is attractive in the scalar and repulsive 
in the pseudoscalar channels, and equal in magnitude. This is the precise behavior needed to lift the degeneracy in 
the constituent models. 

The QCD sum rules offer the possibility to describe mixing and sea quarks effects in glueball states. The method 
is rather technical and requires many ingredients (condensates, instantons, topological charge screening, ...). This 
complexity has led to many different calculations, with different results for the lowest lying states. We refer the reader 
to the literature for the various attempts to describe the low mass glueball spectrum. In here, we present the general 
formalism and some details of a specific calculation, in which the masses are obtained from pure gluonic currents, 
however, sea quarks effects are included (via the quark condensate). 

Glueballs should be produced preferably in gluon- rich processes. We described briefly the main examples: J ftp 
decays, central production pp annihilation,77 fusion and photoproduction. The realization of a new generation of 
experiments, BESIII, PANDA and GlueX, provide hope for new exciting developments in this field. 

The ultra-relativistic heavy ion experiments is another experimental scenario where one expects that glueballs might 
play a role. It is clear that the transition to a different phase of confinement, be it the Strong Coulomb phase or the 
Quark Gluon Plasma, precludes that the confinement properties of the bound systems will change and with it their 
physical behavior. There is no consensus on the glueball properties above the critical temperature and this has led 
to the description of several scenarios which can be very interesting from the point of view of observation. RHIC has 
produced large amount of date at high temperature, as will very soon do ALICE. On the other hand FAIR will probe 
the high density region. Given these circumstances it is clear that the aim is to performed realistic analysis which are 
are able to disentangle the right from the wrong ideas. 

Finally, the strong expected mixing between glueballs and quark states leads to a broadening of the possible glueball 
states which does not simplify their isolation and their theoretical description. The wishful sharp resonances which 
would confer the glueball spectra the beauty, richness and simplicity of the conventional baryonic and mesonic spectra 
are lacking. We hope that in the future and at higher energies the situation changes and we are able to have isolated 
exotic states. It is, however, important to stress, that in any case gluebals are a beautiful and unique consequence of 
QCD. 

The study of glueballs is intimately related to the quantitative understanding of confinement in QCD, since under- 
standing confinement requires an understanding of the soft gluonic field responsible for binding the hadrons and the 
structure of the QCD vacuum. We have emphasized that this has been the goal of the various studies, i.e. to get a 
clearer picture of how QCD behaves at relatively low momentum and how this behavior changes as the temperature 
or the density increase. Our review shows that we have learned much but we have yet not achieved the goal. We 
foresee a bright future in this respect since a new experimental era of the study of QCD in many fronts is opening 
up. To find a description, which is able to describe all the various phenomena in a unique framework, is our task, 
and when achieved will imply that we finally understand QCD. The venture is even more exciting since maybe AdS 
gravity might come to our help, a possibility unthinkable a few years back. 
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